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INTRODUCTION 



In general, random matrices are matrices whose matrix elements are stochastic vari- 
ables. The main goal of Random Matrix Theory (RMT) is to calculate the statistical 
properties of eigenvalues for very large matrices which are important in many appli- 
cations. Ensembles of random matrices first entered in the mathematics literature as a 
p-dimensional generalization of the ^-distribution [1]. Ensembles of real symmetric 
random matrices with independently distributed Gaussian matrix elements were intro- 
duced in the physics literature in order to describe the spacing distribution of nuclear 
levels [2] . The theory of Hermitian random matrices was first worked out in a series of 
seminal papers by Dyson [3]. Since then, RMT has had applications in many different 
branches of physics ranging from sound waves in aluminum blocks to quantum grav- 
ity. For an overview of the early history of RMT we refer to the book by Porter [4]. 
An authoritative source on RMT is the book by Mehta [5]. For a comprehensive review 
including the most recent developments we refer to Ref. [6]. 

Generally speaking, random matrix ensembles provide a statistical description of a 
complex interacting system. Depending on the hermiticity properties of the interactions, 
one can distinguish two essentially different classes of random matrices: Hermitian ma- 
trices with real eigenvalues and matrices without hermiticity properties with eigenvalues 
scattered in the complex plane. We will first give an overview of the ten different classes 
of Hermitian random matrices and then briefly discuss non-Hermitian random matrix 
ensembles. 

The best known random matrix ensembles are the Wigner-Dyson ensembles which are 
ensembles of Hermitian matrices with matrix elements distributed according to 



Here, H is a Hermitian N x N matrix with real, complex, or quaternion real matrix 
elements. The corresponding random matrix ensemble is characterized by the Dyson 
index j8 = 1, 2, and 4, respectively. The measure DH is the Haar measure which is 
given by the product over the independent differentials. The normalization constant of 
the probability distribution is denoted by ,jV . The probability distribution (1) is invariant 
under the transformation 



P(H)DH = J^e 



DH. 



(1) 



H 



UHU 



(2) 
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where U is an orthogonal matrix for /3 = 1, a unitary matrix for /3 = 2, and a symplectic 
matrix for /3 = 4. This is the reason why these ensembles are known as the Gaussian 
Orthogonal Ensemble ( GOE), the Gaussian Unitary Ensemble ( GUE), and the Gaussian 
Symplectic Ensemble ( GSE), respectively. The GOE is also known as the Wishart dis- 
tribution. Since both the eigenvalues of H and the Haar measure DH are invariant with 
respect to (2), the eigenvectors and the eigenvalues are independent with the distribution 
of the eigenvectors given by the the invariant measure of the corresponding orthogonal, 
unitary, or symplectic group. 

There are two ways of arriving at the probability distribution (1). First, from the 
requirement that the matrix elements are independent and are distributed with the same 
average and variance for an ensemble invariant under (2). Second, by requiring that the 
probability distribution maximizes the information entropy subject to the constraint that 
the average and the variance of the matrix elements are fixed. 

A second class of random matrices are the chiral ensembles [7] with the chiral 
symmetries of the QCD Dirac operator. They are defined as the ensembles of N x N 
Hermitian matrices with block structure 

and probability distribution given by 

P(C)DC = ^det N f^ i C o y-^ TlC ' c DC. (4) 

Again, DC is the Haar measure, and Nf is a real parameter (corresponding to the number 
of quark flavors in QCD). The matrix C is a rectangular nx (n + v) matrix. Generically, 
the matrix H in (3) has exactly |v| zero eigenvalues. Also generically, the QCD Dirac 
operator corresponding to a field configuration with the topological charge v has exactly 
| v | zero eigenvalues, in accordance with the Atiyah-Singer index theorem. For this 
reason, v is identified as the topological quantum number. The normalization constant 
of the probability distribution is denoted by JV . Also in this case one can distinguish 
ensembles with real, complex, or quaternion real matrix elements. They are denoted by 
/3 = 1, /3 = 2, and /3 = 4, respectively. The invariance property of the chiral ensembles 
is given by 

C^UCV \ (5) 

where U and V are orthogonal, unitary, and symplectic matrices, respectively. For this 
reason, the corresponding ensembles are known as the chiral Gaussian Orthogonal 
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Ensemble (chGOE), the chiral Gaussian Unitary Ensemble (chGUE), and the chiral 
Gaussian Symplectic Ensemble (chGSE), respectively. A two sub-lattice model with 
diagonal disorder in the chGUE class was first considered in [8]. 

A third class of random matrix ensembles occurs in the description of disordered 
superconductors. Such ensembles with the symmetries of the Bogoliubov-de Gennes 
Hamiltonian have the block structure 



where A is Hermitian and, depending on the underlying symmetries, the matrix B is 
symmetric or anti-symmetric. The probability distribution is given by 



where DH is the Haar measure and jV is a normalization constant. For symmetric B the 
matrix elements of H can be either complex (C) or real (CI). For anti-symmetric B the 
matrix elements of H can be either complex (D) or quaternion real (Dili). The name of 
the ensembles (in parentheses) refers to the symmetric space to which they are tangent 
to. Since they were first introduced by Altland and Zirnbauer [9, 10] we will call them 
the Altland-Zirnbauer ensembles. A hopping model based on the class CI first entered 



A key ingredient in the classification of a Hamiltonian in terms of one of the above 
random matrix ensembles is its anti-unitary symmetries. An anti-unitary operator can 
be written as 



where A is unitary and K is the complex conjugation operator. For the classification 
according to the anti-unitary symmetries we can restrict ourselves to the following three 
different possibilities: (i) the Hamiltonian does not have any anti-unitary symmetries, 
(ii) the Hamiltonian commutes with AK and (AK) 2 = 1, and (iii) [H,AK] = but 
(AK) 2 = -1. In the first case, the matrix elements of the Hamiltonian are complex, 
in the second case, it is always possible to find a basis in which the Hamiltonian is real, 
and in the third case, it can be shown that it is possible to organize the matrix elements 
of the Hamiltonian in quaternion real elements. These three different possibilities are 
denoted by the number of degrees of freedom per matrix element, j8 = 2, j8 = 1, 
and /3 = 4, respectively. This triality characterizes the Wigner-Dyson ensembles, the 
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(6) 




in [11]. 



U=AK, 



(8) 



chiral ensembles, and the Altland-Zirnbauer ensembles. In most cases, the anti-unitary 
operator is the time-reversal symmetry operator. For systems without spin, this is just 
the complex conjugation operator. For systems with spin, the time reversal operator can 

be represented as i(72K, where 02 is one of the Pauli matrices. 

TABLE 1. Random matrix ensemble, corresponding 
symmetric space, and the values for a and j3 . 



RMT 


symmetric space 




a 


GOE 


AI 


1 




GUE 


A 


2 




GSE 


All 


4 




chGOE 


BDI 


1 


\v\+2N f 


chGUE 


ATTT 


2 


l+2|v|+2fy 


chGSE 


CII 


4 


3+4|v|+2A/y 


AZ-CI 


CI 


1 


1 


AZ-D 


D 


2 





AZ-C 


C 


2 


2 


AZ-DIII 


Dili 


4 


1 



We have introduced ten different random matrix ensembles. Each of these ensembles 
can be identified as the tangent space of one of the large families of symmetric spaces 
as classified by Cartan (see Table 1). The matrices in each of these ten ensembles 
can be diagonalized by a unitary transformation, with the unitary matrix distributed 
according to the group measure. For all ensembles, the Jacobian for the transformation 
to eigenvalues as new integration variables depends only on the eigenvalues. For an 
extensive discussion of the calculation of this type of Jacobian we refer to [12]. For the 
Wigner-Dyson ensembles, the joint probability distribution of the eigenvalues is given 
by 

P{{X})d{X} = ^mX})^Y[e- N ^' A dX k , (9) 

where the Vandermonde determinant is defined by 

A({X}) = l\(X k -X l ). (10) 

k>l 

This factor results in correlations of eigenvalues that are characteristic for the random 
matrix ensembles. For example, one finds repulsion of eigenvalues at small distances. 

For the remaining ensembles, the eigenvalues occur in pairs ±\. This results in the 
distribution 

P({X})d{X} = ^\A({X 2 })\Pl\X k a e- N W/ 4 dX k . (11) 
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The values of /3 and a are given in Table 1 below. 

Another well-known random matrix ensemble which is not in the above classification 
is the Poisson Ensemble defined as an ensemble of uncorrelated eigenvalues. Its proper- 
ties are very different from the above RMTs where the diagonalization of the matrices 
leads to strong correlations between the eigenvalues. 

The physical applications of RMT have naturally biased the interest of researchers to 
Hermitian matrices (e.g., the Hamiltonian of a quantum system is a Hermitian operator 
and should be represented by a Hermitian matrix). A variety of methods, described in 
this article, have been developed to treat ensembles of Hermitian matrices. In contrast, 
non-Hermitian random matrices received less attention. Apart from the intrinsic mathe- 
matical interest of such a problem, a number of physically important applications exist 
which warrant the study of non-Hermitian random matrices. 

The simplest three classes of non-Hermitian random matrices, introduced by Ginibre 
[13], are direct generalizations of the GOE, GUE, and GSE. They are given by an 
ensemble of matrices C without any Hermiticity properties and a Gaussian probability 
distribution given by 

P(C)DC = <yfe-^- TrCfc DC , (12) 

where DC is the product of the differentials of the real and imaginary parts of the matrix 
elements of C. Such matrices can be diagonalized by a similarity transformation with 
eigenvalues scattered in the complex plane. The probability distribution is not invariant 
under this transformation, and therefore the eigenvalues and the eigenvectors are not 
distributed independently. Similarly to the Hermitian ensembles, the matrix elements 
can be chosen real, arbitrary complex, or quaternion real. 

The case of the arbitrary complex non-Hermitian random matrix ensemble (12) 
with /3 = 2 is the simplest. The joint probability distribution of eigenvalues {X} = 
{X\ , . . . , Xm} is given by a formula similar to (9): 

P({X})d{X} = ^\A({X})\ 2 l\e- N ^ 2 dx k dy k , (13) 

k 

where x k = ReA^, y k = Im X k . In the quaternion-real case, the joint probability distri- 
bution can also be written explicitly. In the case of real matrices, the joint probability 
distribution is not known in closed analytical form. 

It is also possible to introduce non-Hermitian ensembles with a chiral structure, but 
such ensembles have received very little attention in the literature and will not be 
discussed. What has received a great deal of attention in the literature are non-Hermitian 
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deformations of the Hermitian random matrix ensembles. Among others, they enter in 
the statistical theory of S -matrix fluctuations [14], models of directed quantum chaos 
[15, 16] and in chiral random matrix models at nonzero chemical potential [17]. The 
latter class of ensembles is obtained from (3) and (4) by making the replacement 

C\ ( C+ifl 



r o J " Vr - o (14) 

This chRMT is a model for the QCD partition function at nonzero chemical potential fi 
and will be discussed in more detail below. 

Random Matrix Theory is a theory to describe the correlations of the eigenvalues of a 
differential operator. The correlation functions can be derived from the joint probability 
distribution. The simplest object is the spectral density 

p(X)=^5(X-X k ). (15) 

k 

The average spectral density, denoted by 

R 1 (X) = (p(X)), (16) 

is obtained from the joint probability distribution by integration over all eigenvalues 
except one. The connected two-point correlation function is defined by 

Pc(Ai,A 2 ) = (p(X 1 )p(X 2 )) - (p(A 1 ))(p(A 2 )). (17) 

In RMT it is customary to subtract the diagonal term from the correlation function and 
to introduce the two point correlation function i? 2 (Ai, A 2 ) defined by 

Ri{hM = (p(X l )p(X 2 )) - 8{X X - X 2 )(p{X)) (18) 

and the two-point cluster function 

T 2 {X h X 2 ) = R l (X l )R l (X 2 )-R 2 (X h X 2 ). (19) 

In general, the k— point correlation function can be expressed in terms of the joint 
probability distribution as 

R k (X\,...,X k ) = j^——J dX k+ i---dX N P N (Xi,...,X N ), (20) 

where we have included a combinatorial factor to account for the fact that spectral 
correlation functions do not distinguish the ordering of the eigenvalues. Similarly, one 
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can define higher order connected correlation functions and cluster functions with all 
lower order correlations subtracted out. For details we refer to Mehta's book [5]. 
Instead of the spectral density one often studies the resolvent defined by 



1 _ 1 1 £ 1 



G(z) = -Tr = - V — (21) 

which is related to the spectral density by 

(p(A)) = -lim-Im(G(A +/£)). (22) 

In the analysis of spectra of complex systems and the study of random matrix theories, 
it has been found that the average spectral density is generally not given by the result for 
the Gaussian random matrix ensembles which has a semi-circular shape. What is given 
by RMT are the correlations of the eigenvalues expressed in units of the average level 
spacing. For this reason one introduces the cluster function 

In general, correlations of eigenvalues in units of the average level spacing are called 
microscopic correlations. These are the correlations that can be described by the N — > °° 
limit of RMT. 

The cluster function (23) has universal properties. In the limit N — > °°, it is invariant 
with respect to modifications of the probability distribution of the random matrix en- 
semble. For example, for the GUE and the chGUE it has been shown that replacing the 
Gaussian probability distribution by a distribution given by the exponent of an arbitrary 
even polynomial results in the same microscopic correlation functions [18, 19]. 

For ensembles in which the eigenvalues occur in pairs ±Ajt, an additional important 
correlation function with universal properties is the microscopic spectral density [20] 
defined by 

p s (u) = lim —^j-rlp ( — ^rr J V (24) 

Related to this observable is the distribution of the smallest eigenvalue which was shown 
to be universal as well [21] . For this class of ensembles, the point A = is a special point. 
Therefore, all correlation functions near A = have to be studied separately. However, 
the microscopic correlations of these ensembles in the bulk of the spectrum are the same 
as those of the Wigner-Dyson ensemble with the same value of /3 . 
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There are two different types of applications of RMT. First, as an exact theory of 
spectral correlations of a differential operator. As an important application we mention 
the study of universal properties in transport phenomena in nuclei [14] and disordered 
mesoscopic systems. In particular, the latter topic has received a great deal of attention 
recently (see [6, 42]). This is the original application of RMT. Second, as a schematic 
model for a complex system. One famous example in the second class is the Anderson 
model [22] for Anderson localization. The properties of this model depend in a critical 
way on the spatial dimensionality of the lattice. Other examples that will be discussed 
in more detail below are models for the QCD partition function at nonzero temperature 
and nonzero chemical potential. 

Random Matrix Theory eigenvalue correlations are not found in all systems. Obvi- 
ously, integrable systems, for example a harmonic oscillator, have very different spectral 
properties. Originally, in the application to nuclear levels, it was believed that the com- 
plexity of the system is the main ingredient for the validity of RMT. Much later it was 
realized that the condition for the presence of RMT correlations is that the corresponding 
classical system is completely chaotic. This so called Bohigas-Giannoni-Schmit conjec- 
ture [23] was first shown convincingly for chaotic quantum billiards with two degrees 
of freedom. By now, this conjecture has been checked for many different systems, and 
with some well-understood exceptions, it has been found to be correct. However, a real 
proof is still absent, and it cannot be excluded that additional conditions may be required 
for its validity. In particular, the appearance of collective motion in complex many -body 
systems deserves more attention in this respect. 

In general, the average spectral density is not given by RMT. Therefore, the standard 
procedure is to unfold the spectrum, i.e., to rescale the spacing between the eigenvalues 
according to the local average eigenvalue density. In practice, this unfolding procedure 
is done as follows. Given a sequence of eigenvalues {A^} with average spectral density 
(p(A)), the unfolded sequence is given by 



The underlying assumption is that the average spectral density and the eigenvalue corre- 
lations factorize. The eigenvalue correlations of the unfolded eigenvalues can be investi- 
gated by means of suitable statistics. The best known statistics are the nearest neighbor 
spacing distribution P(S), the number variance T?(r), and the A3 statistic. The number 
variance is defined as the variance of the number of eigenvalues in an interval of length 
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(25) 



r. The A3 statistic is related to the number variance by 

2 f°° 

A 3 (L) = j4j Q (L 3 -2L 2 r + r 3 )L 2 (r)dr. (26) 

In the analysis of spectra it is essential to include only eigenstates with the same 
exact quantum numbers. Spectra with different exact quantum numbers are statistically 
independent. 

The exact analytical expression of the RMT result for the nearest neighbor spacing 
distribution is rather complicated. However, it is well approximated by the Wigner 
surmise which is the spacing distribution for an ensemble of 2 x 2 matrices. It is given 
by 

P(S) =a^e-h s \ (27) 

where the constants ap and bp can be fixed by the conditions that P(S) is normalized 
to unity, and that the average level spacing is one. The level repulsion at short distances 
is characteristic for interacting systems. For uncorrelated eigenvalues one finds P(S) = 
exp(— S). 

Another characteristic feature of RMT spectra is the spectral stiffness. This is ex- 
pressed by the number variance which, asymptotically for large r, is given by 



2 
j8tt 2 



E V) ~ Tr^gr- (28) 



This should be contrasted with the result for uncorrelated eigenvalues given by £ 2 (r) = r. 

In the analysis of spectra one often relies on spectral ergodicity defined as the equiv- 
alence of spectral averaging and ensemble averaging. This method cannot be used for 
the distribution of the smallest eigenvalues, and one necessarily has to rely on ensemble 
averaging. 

Before proceeding to the discussion of mathematical methods of Random Matrix The- 
ory a comment about the notations should be made. Different conventions for normal- 
izing the variance of the probability distribution appear in the literature. This simply 
amounts to a rescaling of the eigenvalues. For example, in the discussion of orthogonal 
polynomials and the Selberg integral below, introduction of rescaled eigenvalues such 
as or XkVN simplifies the expressions. 
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MATHEMATICAL METHODS I: HERMITIAN MATRICES 



Orthogonal polynomials. 

One of the oldest and perhaps most widely used methods in RMT is based on orthogo- 
nal polynomials. A comprehensive presentation of this method is given in Mehta's book 
[5]. Here, we summarize the most important ingredients, concentrating on the GUE for 
mathematical simplicity. 

We have seen in the introductory chapter that the spectral correlation functions can 
be obtained by integrating the joint probability distribution. The mathematical problem 
consists in performing these integrations in the limit N — > °°. It is convenient to rescale Xk 
and introduce Xk = Xk^/N/2. The main point of the orthogonal-polynomial method is the 
observation that the Vandermonde determinant can be rewritten in terms of orthogonal 
polynomials p n (x) by adding to a given row appropriate linear combinations of other 
rows, 

A({*}) = det [(xj) 1 - 1 ] .. = const x det [pi-i( Xj )] .. . (29) 
Including the Gaussian factor in (9), this yields 

N 



_v2 

e 

with functions %(x) satisfying 



A(W) ]J e X 11 = COnSt X det Wi-^ X j)\j=\,...,N ( 3 °) 



J —i 



f°° _ 2 

dx(p m (x)(p n (x) = / dxe x p m {x)p n {x) = 8 nm . (31) 

) J — oo 



In this case, the orthogonal polynomials are essentially the Hermite polynomials, and 
the <p n are the oscillator wave functions, 

e x2 l 2 ( d\ n 2 

%(x)= / w (-l-) e ■ (32) 
y/2 n n\y/% \ dxj 

The integrals in Eq. (20) can now be performed row by row. The fc-point functions are 
then given by determinants of a two-point kernel, 

R k (x h ...,x k ) =det[K N (x i ,x j )] ij=1 j c . (33) 

The kernel K^(x, y) is given by 

N-l 

K N (x,y) = £ (p„(x)(p n (y) (34) 
Random Matrices 1 1 



which can be evaluated using the Christoffel-Darboux formula. In the large-/V limit, the 
spectral density becomes the famous Wigner semicircle, 

Ri(x) = lim K N (x,x) = — \j2N — x 2 (35) 

A?— »°o 7t 

if x 2 < 2N and zero otherwise. The mean level spacing D(x) = l/R\(x) in the bulk of 
the semicircle thus goes like 1/VN. The R k are universal if the spacing \x — y\ is of the 
order of the local mean level spacing, i.e., we require \x — y\ = rD(x) with r of order 
unity. In this limit, we obtain 

sin i Ttv ) 

D(x) lim K N (x,y) = — — (36) 

which is the famous sine kernel. The various functions appearing in a typical RMT- 
analysis, e.g., P(s), £ 2 («), or A^n), can all be expressed in terms of the R^. 



Selberg's integral. 

In 1944, Selberg computed an integral which turned out to have significant applica- 
tions in RMT [24]. His result reads [5] 

I(a,p,y,n) = dx, ■ ■ ■ j\x n \A{x)\ 2 vf\x<*-\\ -xj)^ 1 

fr 1 r(l + y+jy)r(« + j7)r(/3+jy) 
" fj r(l + y)r(a + p + (n + j-l)y) ' 1 ; 

where A(x) is the Vandermonde determinant, n is an integer, and a, /3, and y are complex 
numbers satisfying Re a > 0, Re /3 > 0, Rey> -min{l/n,Rea/(n - 1), Re/3/(n- 
1)}. Choosing the parameters in Eq. (37) appropriately, one can derive special forms of 
Selberg's integral related to specific orthogonal polynomials [5]. For example, choosing 
Xi =yi/L, a = j8 = aL 2 + 1, and taking the limit L — > one obtains the integrals of the 
joint probability density function of the GUE which are related to Hermite polynomials. 
Selberg's integral is also very useful in the derivation of spectral sum rules [25]. 
Aomoto derived the following generalization of Selberg's integral [26], 

j^dxx ■■■ j\x n xv--x m \A{x)\ 2 vf\x c ;-\\-x j f-' 
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where 1 < m < n. A further extension of Selberg's integral was considered by Kaneko 
[27] who related it to a system of partial differential equations whose solution can be 
given in terms of Jack polynomials. 



Supersymmetric method 

The supersymmetric method has been applied successfully to problems where the 
orthogonal polynomial method has failed [28, 29, 14]. It relies on the observation that 
the average resolvent can be written as 



,W XX 1 Irr, 1 \ 



N\ z-H I N dJ 
where the generating function is defined by 



Z(J), (39) 



and the integral is over the probability distribution of one of the random matrix ensem- 
bles defined in the introductory chapter. The determinant can be expressed in terms of 
Gaussian integrals, 



J d{ ¥ } e xp(^[i^(z-H) kl <j> l + ixt(z + J-H) kl Xi}^ , (41) 



det(z-H + J) 
det(z-H) 

where the measure is defined by 

w-n'ww . (42) 

For convergence the imaginary part of z has to be positive. The integrations over the real 
and imaginary parts of fa range over the real axis (the usual commuting, or bosonic vari- 
ables), whereas Xi an d Xi are Grassmann variables (i.e., anticommuting, or fermionic 
variables) with integration defined according to the convention that 



dx = and / xdX = 1 • (43) 



With this normalization, Z(0) = 1 . 

For simplicity, we consider the GUE [/3 =2 in (1)] which mathematically is the 
simplest ensemble. The Gaussian integrals over H can be performed trivially, resulting 
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in the the generating function 



Z(J) 



Jd{ ¥ } 



exp 



Trg 

IN 5 



i 



(44) 



where the sums over j run from 1 to N. The symbol Trg denotes the graded trace (or 
supertrace) defined as the difference of the trace of the boson-boson block (upper left) 
and the trace of the fermion-fermion (lower right) block. For example, in terms of the 
2x2 matrix (46), Trgc = Obb — i<Jff- The quartic terms in and % can be expressed as 
Gaussian integrals by means of a Hubbard-Stratonovitch transformation. This results in 



Z(J) = J d{y}doexp 



N 



where 



--Tr gC 7 2 + /£ 



(45) 



and 



z 
z + J 



(46) 



(47) 



The variables Obb and Off are commuting (bosonic) variables that range over the full 
real axis. Both Obf and Ofb are Grassmann (fermionic) variables. 

The integrals over the <j> and the % variables are now Gaussian and can be performed 
trivially. This results in the a-model 



Z(7) = J do exp 



■^Trg(7 2 +iVTrglog((7 + C) 



By shifting the integration variables according to o 
differentiation with respect to J one easily finds that 

(G(z)) = (z-ia FF ). 



(48) 

o — £ and carrying out the 

(49) 



In the large N limit, the expectation value of Off follows from a saddle-point analysis. 
The saddle point equation for Off is given by 
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resulting in the resolvent 

{G(z)) = Z -- l -VA^. (51) 

Using the relation (22) we find that the average spectral density is a semi-circle. 

The super symmetric method can also be used to calculate spectral correlation func- 
tions. They follow from the average of the advanced and the retarded resolvent. In that 
case we do not have a saddle-point but rather a saddle-point manifold related to the hy- 
perbolic symmetry of the retarded and advanced parts of the generating function. The 
supersymmetric method not only provides us with alternative derivations of known re- 
sults. As an example we mention that the analytical result for S-matrix fluctuations at 
different energies was first derived by means of this method [14]. 

Alternatively, it is possible to perform the o integrations by a supersymmetric version 
of the Itzykson-Zuber integral [30] rather than a saddle-point approximation. The final 
result is an exact expression for the kernel of the correlation functions. The advantage 
of this method is that it exploits the determinantal structure of the correlation functions 
[see (33)] and all correlations functions are obtained at the same time. Moreover, the 
results are exact at finite N. 



Replica trick 

The replica trick, which was first introduced in the theory of spin glasses [31], is 
based on the observation that 



v 1 \ v 1 d 

Tr ) = hm - — 

z — Hf n^on dJ 



Zn(J), (52) 

7=0 



where the generating function is defined by 

Z n (J) = J DHP(H)det n (z-H + J). (53) 

The determinant can then be expressed as a Grassmann integral where the ^-variables 
now have an additional flavor index, 

dot n (z-H + J) = [ d{ X }expY,x(*(z-H + J) ijX f. (54) 
The sum over / ranges from 1 to n, and the measure is defined by 

d{x}=f\\\dx{dxt- (55) 

f=\i=\ 
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After averaging over the matrix elements of H and a Hubbard-Stratonovitch transforma- 
tion one can again proceed to the o variables. In this case, we only have a Off block 
which is now an n x n matrix. The average resolvent then follows by making a saddle 
point approximation and taking the replica limit with the same final result as given in 
(51). 

Because the replica trick relies on an analytical continuation in n, it is not guaranteed 
to work. Several explicit examples for its failure have been constructed [32, 17]. In 
general, it cannot be used to obtain nonperturbative results for eigenvalue correlations 
on the microscopic scale which decreases as l/N in the limit N — > °°. 



Resolvent Expansion Methods 

The Gaussian averages can also be performed easily by expanding the resolvent in a 
geometric series in l/z, 

G(z) = - + i^r-H- + Ltv-H-H- + ■■■. (56) 

zNzzNzzz 

The Gaussian integral over the probability distribution of the matrix elements is given 
by the sum over all pairwise contractions. For the GUE, a contraction is defined as 

(HijHl) = ^8 a 8 jk . (57) 

To leading order in l/N, the contributions are given by the nested contractions. One 
easily derives that the average resolvent satisfies the equation 

(G(z)) = - z (\ + (G(z)) 2 ), (58) 

again resulting in the same expression for the average resolvent. 

This method is only valid if the geometric series is convergent. For this reason, the 
final result is only valid for the region that can be reached from large values of z by 
analytical continuation. For non-Hermitian matrices this leads to the failure of this 
method, and instead one has to rely on the so-called Hermitization. 

As is the case with the replica trick, this method does not work to obtain nonperturba- 
tive results for microscopic spectral correlations. This method has been used widely in 
the literature. As one of the earlier references we mention the application to the statistical 
theory of nuclear reactions [33] 
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Dyson Gas 



The formula (9) suggests a very powerful analogy between the Wigner-Dyson random 
matrix ensembles and the statistical properties of a gas of charged particles restricted to 
move in one dimension, the Dyson gas [3]. 

Let Xk be a coordinate of a classical particle which moves in the potential Vi(A^) = 
NX?/4. Furthermore, let two such particles repel each other so that the potential of 
the pairwise interaction is V^(Ajt,A;) = — ln|A^ — A/|. If one considers a gas of N 
such particles in thermal equilibrium at temperature T, then the probability distri- 
bution for the coordinates of the particles A = will be proportional to 
exp(— V (A)/r) Ylkd^k, where the potential energy V is given by 

^(^) = E^2(A fc ,A / )+£yi(A ; -). (59) 

k<l i 

If the temperature T of the gas is chosen to be equal to 1//3, the probability distribution 
of the coordinates of the particles becomes identical to the probability distribution (9) of 
the eigenvalues. 

This analogy allows one to apply methods of statistical mechanics to calculate dis- 
tributions and correlations of the eigenvalues [3]. It also helps to grasp certain aspects 
of universality in the statistical properties of the eigenvalues. In particular, it is under- 
standable that the correlations in the relative positions of particles are determined by the 
interactions between them, i.e., by V%, and are generally insensitive to the form of the 
single-particle potential, V\ . On the other hand, the overall density will depend on the 
form of the potential V\ . 

The logarithmic potential V2 is the Coulomb potential in two-dimensional space (i.e., 
it satisfies the two-dimensional Laplace equation AV2 = 0). Therefore, the Dyson gas 
can be viewed as a two-dimensional Coulomb gas, with the kinematic restriction that 
the particles move along a straight line only. This restriction is absent in the case of 
non-Hermitian matrices. 



MATHEMATICAL METHODS II: NON-HERMITIAN MATRICES 

The eigenvalues of non-Hermitian matrices are not constrained to lie on the real axis. 
Instead, they occupy the two-dimensional complex plane. This fact requires non-trivial 
modifications of some of the methods developed for Hermitian matrices. Surprisingly or 
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not, the required formalism is sometimes simpler, and sheds more light on the properties 
of Hermitian random matrices. 



Orthogonal polynomials 

The method of orthogonal polynomials can also be applied to treat non-Hermitian 
random matrices. The simplest example is the Ginibre ensemble of arbitrary complex 
matrices (12) with /3 = 2 [13]. It is convenient to rescale X k and introduce w k = A^a/N. 
The orthogonal polynomials with respect to the weight given by exp(— \w\ 2 ) are simply 
the monomials w n . Indeed, 

J dudve- M2 w n (w*) m = Kn\8 mn , (60) 

where w = u + iv. The orthonormal functions, / dudv(j) n (w)<j> m (w*) = 5 nm , are, therefore, 

<j> n (w) = — Le-M 2 ^. (61) 
V Kn\ 

Following the same steps as in the case of the Hermitian GUE one obtains all correlation 
functions in the form of the determinant 

R k (w h ...,w k ) = det[K N (wi,Wj)] iJ= i^ k , (62) 

with a kernel Kpj given by 

N-l 

K N (w h W 2 ) = <l>n{w\)<l>m{wl). (63) 



By a careful analysis of the large-Af limit of the kernel one finds that7?i (w) is \/iz inside 
the complex disk \w\ < VN and vanishes outside this domain. 



Coulomb gas 

The probability distribution (13) is the same as for a Coulomb gas in two dimensions 
placed in the harmonic potential V\ = N\z\ 2 /2 = N(x 2 +y 2 )/2 at a temperature 1//3 = 
1/2. Unlike in the Hermitian case, the particles of the gas are now allowed to move in 
both dimensions. 

The analogy with the Coulomb gas can be used to calculate the density of eigenvalues 
of the ensemble of complex non-Hermitian matrices [(12) with /3 = 2] in the limit 
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N — > oo. In this limit, the typical energy per particle, 0(N), is infinitely larger than 
the temperature, 1/2. Therefore, the system is assuming an equilibrium configuration 
with the minimal energy, as it would at zero temperature. Each particle is subject to 
a linear force —dV\/d\z\ = —N\z\ directed to the origin, z = 0. This force has to be 
balanced by the Coulomb forces created by the distribution of all other particles. Thus, 
the electric field created by this distribution must be directed along the radius and be 
equal to \E\ = N\z\. The Gauss law, VE = 2np, tells us that such a field is created 
by charges distributed uniformly (with density p = N/k) inside a circle around z = 0, 
known as the Ginibre circle. The radius of this circle, R, is fixed by the total number of 
the particles, 7tR z p = N, so that R= 1. 



Electrostatic analogy and analyticity of resolvent 

In general, the mapping of the random matrix model onto the Coulomb gas is not 
possible, because the pairwise interaction is not always given simply by the logarithm 
of the distance between the particles. However, a more generic electrostatic analogy 
exists, relating the two-dimensional density of eigenvalues p, 

p(x,y) = Y, 8 ( x - x k) 8 (y-yk), (64) 

k 

where x^ and y k are real and imaginary parts of X k , and the resolvent G, 

Since the electric field created by a point charge in two dimensions is inversely pro- 
portional to the distance from the charge, one can see that the two-component field 
(iVReG, — MmG) coincides with the electric field E, created by the charges located at 
the points {Ai , . . . , Ayy} in the complex plane. 

The Gauss law, relating the density of charges and the resulting electric field, VE — 
2k p, gives the following relation between the density of the eigenvalues and the resol- 
vent: 

_ N ( dReG dlmGl 
^ In \ dx dy J 
This relation is the basis of methods for the calculation of the average density of the 
eigenvalues, (p). The r.h.s of this equation vanishes if G obeys the Cauchy-Riemann 
conditions, i.e., if it is an analytic function of the complex variable z = x + iy. Conversely, 
p describes the location and the amount of non-analyticity in G. 

Random Matrices 19 



In the case of Hermitian matrices, C = H, the eigenvalues lie on a line (real axis) and, 
after ensemble averaging, they fill a continuous interval. This means that the average 
resolvent (G(x,y)) has a cut along this interval on the real axis. The discontinuity along 
this cut is related to the linear density of the eigenvalues by (22). In the case of a 
non-Hermitian matrix C, the eigenvalues may and, in general, do fill two-dimensional 
regions. In this case, the function G is not analytic in such regions. 

This is best illustrated by the ensemble of arbitrary complex matrices C (12). In the 
N — > °° limit, the resolvent is given by 

= { f/'z, N > ! ■ <67) 

One observes that G is non-analytic inside the Ginibre circle [13]. 



Replica trick 

The generalization of the replica trick to the case of non-Hermitian matrices is based 
on the relation 

/tt-!-\ = lim-|-lnZ M (z), (68) 
\ z — C I n^O n dz 

where now 

Z n (z) = J DCP(C)\det(z-C)\ n . (69) 

The absolute value of the determinant can be also written as det"/ 2 (z - C) deW 2 (z* - C t ) . 
Following (54), one introduces n/2 Grassmann variables Xi to represent deW 2 (z — C) 
and another n/2 to represent det w / 2 (z* — C"). If the measure P(C) is Gaussian, the 
integral over C can now be performed resulting in terms quartic in the Grassmann 
variables. These can be rewritten with the help of an auxiliary n x n variable o as 
bilinears in %, after which the % integration can be done. The resulting integral, in the 
limit N — > oo, is given by the saddle point (maximum) of its integrand. In the case of the 
Ginibre ensembles one arrives at the following expression: 

lnZ n (z) =^max[-|(7| 2 + ln(|z| 2 + |c7| 2 )l . (70) 
a 

There are two possible maxima, a = and \o\ 2 = 1 — |z| 2 , which give two branches for 
lnZ M , lnZ„/ (nN) = In \z\ 2 and lnZ„/ (nN) = \z\ 2 — 1. The former dominates when |z| > 1 
and the latter when |z| < 1. Using (68) one obtains the average resolvent given by (67). 
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It is important that the absolute value of the determinant is taken in (69). Without 
taking the absolute value one would obtain the incorrect result G= l/z everywhere in 
the complex plane. 

Hermitization 

The method of Hermitization, as well as the replica trick (68,69), is based on the 
observation that Aln |z| 2 = 4n8(x)8(y), where A is the Laplacian in the coordinates x 
and y. One can, therefore, write for the eigenvalue density p(x,y) 

4n(p (x,y) ) = A(ln det(z - C) {z* - C f ) ) . (71) 

The determinant on the r.h.s. can be written as the determinant of a matrix (up to a sign) 

tf(z,z*)=( z *.° c t Z ~°y (72) 

This matrix is Hermitian, and one can apply methods of Hermitian RMT (e.g., the 
supersymmetric method or the replica trick) to determine its resolvent G(rj). Integrating 
over r\ one obtains the quantity 

(lndet(T]-tf)), (73) 
which in the limit r\ — > reduces to the expression on the r.h.s. of (71) [34, 15, 35, 36]. 

APPLICATIONS AND ADVANCED TOPICS 

In this section, we briefly review a variety of different subfields of physics where RMT 
has been applied successfully. Most of the examples can be found in the comprehensive 
presentation of Ref . [6] which also contains a wealth of useful references. 

Nuclear level spacings. 

Historically, the first application of RMT in physics arose in the study of nuclear 
energy levels. The problem of computing highly excited energy levels of large nuclei is 
so complicated that it is impossible to make detailed predictions based on microscopic 
models. Therefore, as discussed in the introduction, it is interesting to ask whether 
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FIGURE 1. The histogram represents the nearest-neighbor spacing distribution of the "Nuclear Data 
Ensemble" (NDE). The curve labeled GOE is the random-matrix prediction, and the Poisson distribu- 
tion, representing uncorrected eigenvalues, is shown for comparison. Taken from Ref. [37] with kind 
permission from Kluwer Academic Publishers. 

the statistical fluctuations of the nuclear energy levels are universal and described by 
the predictions from RMT. The nuclear Hamiltonian is time-reversal invariant so that 
the data should be compared with GOE results. Figure 1 shows the nearest neighbor 
spacing distribution of nuclear energy levels of the same spin and parity, averaged over 
1726 spacings from 32 different nuclei [37]. Clearly, the data are described by RMT, 
indicating that the energy levels are strongly correlated. The parameter-free agreement 
seen in the figure gave strong support to the ideas underlying RMT. 

Hydrogen atom in a magnetic field. 

The Hamiltonian of this system is given by 

H = |— - — — — G)L Z + \co\x 2 + y 2 ) , (74) 
2m r 2 

where m is the reduced mass, e is the unit charge, r = (x 2 + y 2 +z 2 ) 1//2 is the separation 
of proton and electron, co = eB/ (2mc) is the Larmor frequency, B is a constant magnetic 
field in the z-direction, and L z is the third component of the angular momentum. At 
5 = 0, the system is integrable. This property is lost when the magnetic field is turned on, 
and large parts of the classical phase space become chaotic. For an efficient numerical 
computation of the eigenvalues, it was important to realize that the Hamiltonian has a 
scaling property which simplifies the calculations considerably: The spectrum depends 
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only on the combination £ = y~ 2 l 3 E, where 7 is a dimensionless variable proportional 
to B and E is the energy of the system [38]. This variable increases if the magnetic 
field is increased and/or the ionization threshold is approached from below. Thus, as 
a function of £, one should observe a transition from Poisson to RMT behavior in the 
spectral correlations. The numerical results are in agreement with experimental data and 
clearly show a Poisson to RMT transition, see Fig. 2. 
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FIGURE 2. Nearest neighbor spacing distribution of the energy levels of the hydrogen atom in a 
magnetic field (histograms). The solid line in the bottom plot is the RMT prediction for the GOE, all 
other lines are fits. As a function of the scaled variable e, which increases from top to bottom, a transition 
from Poisson [P(x) = exp(— x)] to RMT behavior is observed. Taken from Ref. [39] with kind permission 
from Elsevier Science. 



Billiards and quantum chaos. 



These are the prototypical systems used in the study of quantum chaos. A billiard 
is a dynamical system consisting of a point particle which can move around freely in a 
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bounded region of space. For simplicity, we assume that the space is two-dimensional. In 
a classical billiard, the particle is reflected elastically from the boundaries corresponding 
to a potential that is zero inside the boundary and infinite outside the boundary. In a 
quantum billiard this results in a free particle Schrddinger equation with wave functions 
that vanish on and outside this boundary. Depending on the shape of the boundary, 
the classical motion of the particle can be regular or chaotic (or mixed). Examples 
of classically regular billiards are the rectangle and the ellipse. Important classically 
chaotic billiards are the stadium billiard (i.e., two semicircles at opposite sides of an open 
rectangle) and the Sinai billiard (i.e., the region outside a circle but inside a concentric 
square surrounding the circle). According to the conjecture by Bohigas, Giannoni, and 
Schmit [23], the level correlations of a quantum billiard whose classical counterpart is 
chaotic should be given by RMT, whereas the eigenvalues of a quantum billiard whose 
classical analog is regular should be uncorrelated and, thus, described by a Poisson 
distribution. This conjecture was investigated — numerically, semiclassically, or using 
periodic orbit theory — in a number of works and confirmed in almost all cases [40]. One 
can also vary the shape of a billiard as a function of some parameter, thus interpolating 
between a classically regular and a classically chaotic billiard. As a function of the 
parameter, one then observes a transition from Poisson to RMT behavior in the level 
correlations of the corresponding quantum billiard. 



Quantum dots. 

Semiconducting micro structures can be fabricated such that the electrons are confined 
to a two-dimensional area. If this region is coupled to external leads, we speak of a 
quantum dot. Such systems have many interesting properties. If the elastic mean free 
path of the electrons (which at very low temperatures is ^10/im) is larger than the linear 
dimensions (~ 1 /im) of the quantum dot, and if the Coulomb interaction is neglected, 
the electrons can move around freely inside the boundary, and the quantum dot can be 
thought of as a realization of a quantum billiard. Depending on the shape of the quantum 
dot, certain observables, e.g., the fluctuations of the conductance as a function of an 
external magnetic field, show a qualitatively different behavior. If the shape is classically 
chaotic (e.g., a stadium), the experimental results agree with predictions from RMT as 
expected, in contrast to data obtained with quantum dots of regular shape where the 
fluctuations are not universal [41]. For a recent review of quantum dots and universal 
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conductance fluctuations to be discussed in the following section we refer to Ref. [42]. 

Universal conductance fluctuations. 

A mesoscopic system in condensed matter physics is a system whose linear size 
is larger than the elastic mean free path of the electrons but smaller than the phase 
coherence length, which is essentially the inelastic mean free path. A typical size is of 
the order of 1 /im. The conductance, g, of mesoscopic samples is closely related to their 
spectral properties. Using a scaling block picture, Thouless found that in the diffusive 
regime, g = Ec/A, where Ec/his the inverse diffusion time of an electron through the 
sample and A is the mean level spacing [43] . This can be rewritten as g = (N(Ec) ) , where 
(N(E)) is the mean level number in an energy interval E. Thus the variance, (8g 2 ), of 
the conductance is linked to the number variance, E 2 , of the energy levels. 

In experiments at very low temperatures where the conductance of mesoscopic wires 
was measured as a function of an external magnetic field, people have observed fluctua- 
tions in g of the order of e 2 /h, independent of the details of the system (shape, material, 
etc.). These are the so-called universal conductance fluctuations [44]. This phenomenon 
can be understood qualitatively by estimating the number fluctuations of the electron 
levels using RMT results. However, the magnitude of the effect is much larger than 
expected, due to complicated quantum interference effects. While a truly quantitative 
analysis requires linear response theory (the Kubo formula) or the multichannel Lan- 
dauer formula, both the magnitude of the fluctuations as well their universality can be 
obtained in a simpler approach using the transfer matrix method. Here, the assumption 
(although not quantitatively correct) is that certain parameters of the transfer matrix have 
the same long-range stiffness as in RMT spectra. 

Anderson Localization 

Anderson localization is the phenomenon that a good conductor becomes an insulator 
when the disorder becomes sufficiently strong. Instead of a description of the electron 
wave functions by Bloch waves, the wave function of an electron becomes localized and 
decays exponentially, i.e., 

\jf(r) ~e- r/Lc . (75) 
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The length scale L c is known as the localization length. This phenomenon was first 
described in the Anderson model [22] which is a hopping model with a random potential 
on each lattice point. The dimensionality of the lattice plays an important role. It has 
been shown that in one dimension all states are localized. The critical dimension is two, 
whereas for d = 3 we have a delocalization transition at an energy El- All states below 
El are localized whereas all states above El are extended, i.e., with a wave function that 
scales with the size of the system. 

The eigenvalues of the localized states are not correlated, and their correlations are 
described by the Poisson distribution. In the extended domain the situation is more 
complicated. An important energy scale is the Thouless energy [43] which is related 
to the diffusion time of an electron through the sample. With the latter given by L 2 /D 
(the diffusion constant is denoted by D) this results in a Thouless energy given by 

hD 

E c = jj. (76) 

Correlations on an energy scale below the Thouless energy are given by Random Matrix 
Theory, whereas on higher energy scales the eigenvalues show weaker correlations. 



Other wave equations. 

So far, we have implicitly considered quantum systems which are governed by the 
Schrodinger equation. It is an interesting question to ask if the eigenmodes of systems 
obeying classical wave equations display the same spectral fluctuation properties as 
predicted by RMT Classical wave equations arise, e.g., in the study of microwave 
cavities or in elastomechanics and acoustics. 

In three-dimensional microwave cavities, the electric and magnetic fields are deter- 
mined by the Helmholtz equation, (A 2 + k 2 )A(r) = 0, where A = E or B. It was found 
experimentally that the spacing of the eigenmodes of the system is of RMT type if the 
cavity has an irregular shape [45]. If the cavity has some regular features, the spacing 
distribution interpolates between RMT and Poisson behavior [46]. 

Elastomechanical eigenmodes have been studied both for aluminum and for quartz 
blocks. Here, there are two separate Helmholtz equations for the longitudinal (pressure) 
and transverse (shear) waves, respectively, making the problem even more different 
from the Schrodinger equation. Several hundred up to about 1500 eigenmodes could 
be measured experimentally. A rectangular block has a number of global symmetries, 
and the measured spectrum is a superposition of subspectra belonging to different 
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symmetries. In such a situation, the spacing distribution of the eigenmodes is expected 
to be of Poisson type, and this was indeed observed experimentally. The symmetry can 
be broken by cutting off corners of the block, and the resulting shape is essentially a 
three-dimensional Sinai billiard. Depending on how much material was removed from 
the corners, a Poisson to RMT transition was observed in the spacing distribution of the 
eigenmodes [47]. 

Thus, we conclude that RMT governs not only the eigenvalue correlations of the 
Schrodinger equation but also those of rather different wave equations. 



This is an example from number theory which, at first sight, is not related to the 
theory of dynamical systems. The Riemann zeta function is defined by £(z) = Y*k=\ ^ Z 
for Rez > 1. Its nontrivial zeros z n are conjectured to have a real part of 1/2, i.e., 
z n = 1/2 + iy n - An interesting question is how the y n are distributed on the real axis. 
To this end, it was argued that the two-point correlation function of the y n has the form 
Yi{ r ) = 1 — [sin(TTr) / (nr)] 2 [48]. This is identical to the result obtained for the unitary 
ensemble of RMT and consistent with a conjecture (apparently by Polya and Hilbert) 
according to which the the zeros of £(z) are related to the eigenvalues of a complex 
Hermitian operator. By computing the y n numerically up to order 10 20 [49], it was shown 
that their distribution indeed follows the RMT prediction for the unitary ensemble (for 
large enough y„). 

Universal eigenvalue fluctuations in Quantum chromodynamics 



QCD is the theory of the strong interactions, describing the interaction of quarks 
and gluons which are the basic constituents of hadrons. QCD is a highly complex and 
nonlinear theory for which most nonperturbative results have been obtained numerically 
in lattice QCD using the worlds fastest supercomputers. The Euclidean QCD partition 
function is given by 



where Sym is the Euclidean Yang-Mills action and the path integral is over all SU (N c ) 
valued Yang-Mills fields A l „ (pi is the Lorentz index, N c the number of colors, and 



Zeros of the Riemann zeta function. 



(QCD) 




(77) 
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Nf the number of quark flavors). The Euclidean Dirac operator is defined by D = 
Yfidfi + igYfiA^, where g is the coupling constant and 7^ are the Euclidean gamma 
matrices. Because of the chiral symmetry of QCD, in a chiral basis the matrix of D 
has the block structure 

a>=(£ (78) 

In a lattice formulation, the dimension of the matrix T is a multiple of the total number 
of lattice points. The smallest eigenvalues of the Dirac operator play an important role 
in the QCD partition function. In particular, the order parameter of the chiral phase 
transition is given by 

E=limlim 5<£M, (79) 

where (p(A)) is the average spectral density of the Dirac operator and V is the volume 
of space-time. 

Although the QCD partition function can only be calculated numerically, in certain 
domains of the parameter space it is possible to construct effective theories which can 
be solved analytically. An important ingredient is the chiral symmetry of the QCD 
Lagrangian which is broken spontaneously in the ground state. Considering Euclidean 
QCD in a finite volume, the low-energy behavior of the theory can be described in terms 
of an effective chiral Lagrangian if the linear length, L, of the box is much larger than 
the inverse of a typical hadronic scale. Furthermore, if L is smaller than the inverse of 
the mass of the pion, which is the Goldstone boson of chiral symmetry breaking, then 
the kinetic terms in the chiral Lagrangian can be neglected. It was found by Leutwyler 
and Smilga that the existence of this effective partition function imposes constraints on 
the eigenvalues of the QCD Dirac operator [50]. However, in order to derive the full 
spectrum of the Dirac operator one needs a different effective theory defined by the 
partially quenched chiral Lagrangian, which in addition to the usual quarks includes a 
valence quark and its superpartner [51]. As is the case with the usual chiral Lagrangian, 
the kinetic terms of this Lagrangian can be neglected if the inverse mass of the Goldstone 
bosons corresponding to the valence quark mass is much larger than the size of the box. 
It has been shown that in this domain the corresponding spectral correlators are given 
by the chiral ensembles which have the same block structure as the Dirac operator (78). 
The j8 -value of the ensemble is determined by N c and the representation of the fermions. 
The energy scale for the validity of chiral RMT is the equivalent of the Thouless energy 
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and is given by F 2 /ZL 2 , where F is the pion decay constant that enters in the chiral 
Lagrangian. 

The fluctuation properties of the Dirac eigenvalues can be studied directly by diago- 
nalizing the lattice QCD Dirac operator. Correlations in the bulk of the spectrum agree 
perfectly with the various RMT results [52]. However, as was already pointed out, the 
small Dirac eigenvalues are physically much more interesting. Because of the relation 
(79), the spacing of the low-lying eigenvalues goes like 1/(VX). To resolve individual 
eigenvalues one has to magnify the energy scale by a factor of VE and consider the 
microscopic spectral density [20] defined in (24). 

Because of the chiral structure of the Dirac operator in (78), all nonzero eigenvalues 
of iD come in pairs ±A„, leading to level repulsion at zero. This is reflected in the fact 
that p s (0) = even though \imi^ t0 \im v ^ oo p(?i) /V > 0. The spectrum is said to have a 
"hard edge" at A = 0. 
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FIGURE 3. Distribution of the smallest eigenvalue (left) and microscopic spectral density (right) of the 
QCD Dirac operator. The histograms represent lattice data in quenched SU(2) with staggered fermions 
on a 10 4 lattice using j3 = A/g 2 = 2.0 (not to be confused with the Dyson index j3). The dashed curves 
are the parameter-free RMT predictions. Taken from Ref. [54] with kind permission from the American 
Physical Society. 

The result for p s (z) for the chGUE (appropriate for QCD with three and more colors) 
and gauge fields with topological charge v reads [53, 7] 



2 

Ps{z) = -z J Nf+ \ v \(z) ~~ ^V/+|v| + l( z )^V/+|v|-l(^) 



(80) 



2 Y N s A 

where J denotes the Bessel function. The results for the chGSE and the chGOE are more 
complicated. Lattice QCD data agree with RMT predictions as seen in Fig. 3 which 
represents results corresponding to the chGSE. 
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QCD at nonzero temperature and chemical potential 



Random-matrix models can also be used to model and analyze generic properties of 
the chiral symmetry restoration phase transition at finite temperature or finite baryon 
chemical potential /i. For example, the effect of the chemical potential can be described 
by the non-Hermitian deformation (14) of the chGUE. The eigenvalues of such a matrix 
are not constrained to lie on the real axis. The quantity that signals chiral symmetry 
breaking is the discontinuity (a cut) of the averaged resolvent (G(z)) at z = 0. One 
can calculate (G(z)) in a theory with n ^0, which corresponds to QCD with n species 
of quarks. There is a critical value of /i above which (G(z)) becomes continuous at 
z = 0, and, therefore, chiral symmetry is restored. In lattice Monte Carlo the problem 
of calculating the partition function and expectation values such as (G(z)) at finite /i, 
which are of a paramount interest to experiment, is still unresolved. The difficulty lies in 
the fact that the determinant of the Dirac matrix is complex and cannot be used as part 
of the probabilistic measure to generate configurations using the Monte Carlo method. 
For this reason, exploratory simulations at finite /i have only been done in the quenched 
approximation in which the fermion determinant is ignored, n = 0. The results of such 
simulations were in a puzzling contradiction with physical expectations: The transition 
to restoration of chiral symmetry occurs at /i = in the quenched approximation. 

The chiral random matrix model at ji ^ allows for a clean analytical explanation 
of this behavior, since one can calculate (G(z)) both at n = and n ^ 0. (As before, 
the number of replicas is denoted by n.) The behavior of (G(z)) at n = and n ^ 
is drastically different. While at n ^ the non-analyticities of (G(z)) come in the 
form of one-dimensional cuts, for n = they form two-dimensional regions, similar 
to the Ginibre circle in the case of the non-Hermitian GUE. This means that the n = 
(quenched) theory is not a good approximation to the n ^ (full) theory at finite /i, 
when the Dirac operator is non-Hermitian. The quenched theory is an approximation 
(or, the n — > limit) of a theory with the determinant of the Dirac operator replaced by 
its absolute value, which has different properties at finite fi [17]. 



Quantum gravity in two dimensions. 

In all cases we have discussed so far, the random-matrix model was constructed for 
the Hamiltonian (or a similar operator) of the system, and the universal properties were 
independent of the distribution of the random matrix. In contrast, in quantum gravity 



Random Matrices 



30 



the elementary fields are replaced by matrices, and the details of the matrix potential do 
influence the results. For a recent review we refer to [55]. 

Two dimensional quantum gravity is closely related to string theory. The elementary 
degrees of freedom are the positions of the string in d dimensions. The action, S, of the 
theory involves kinematic terms and the metric. The partition function, Z, is then given 
as a path integral of exp(— S) over all possible positions and metrics. The string sweeps 
out two-dimensional surfaces, and Z can be computed in a so-called genus expansion, 
i.e., as a sum over all possible topologies of these surfaces. This is typically done by 
discretizing the surfaces. One can then construct dual graphs by connecting the centers 
of adjacent polygons (with n sides). These dual graphs turn out to be the Feynman 
diagrams of a <p n -theory in zero dimensions which can be reformulated in terms of a 
matrix model. The partition function of this model is given by 



where the M are Hermitian matrices of dimension N and the g n are coupling constants 
involving appropriate powers of the cosmological constant. The mathematical methods 
used to deal with the matrix model of quantum gravity are closely related to those 
employed in RMT, giving rise to a useful interchange between the two areas. 
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